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ABSTRACT 


This program alms at developing mathematical models and computer codes 
based on these models, which allow prediction of the product distribution In 
chemical reactors for converting gaseous silicon compounds to condensed-phase 
silicon. The major Interest Is In collecting silicon as a liquid on the reac- 
tor walls and other collection surfaces. Two reactor systems are of major 
Interest, the Uestlnghouse SlCl<./Na reactor In which Si(0 is collected on the 
flow tube reactor walls and the new AeroChem reactor in which Si(l) droplets 
formed by the SiClt/Na reaction are collected by a jet Impingement method. 

During this quarter the following tasks have been accomplished: Cl) 

formulation of a model for silicon vapor separation/collection from the de- 
veloping turbulent flow stream within reactors of the Westinghouse type, (11) 
modification of an available general parabolic code to achieve solutions to 
the governing partial differential equations (boundary layer type) which de- 
scribe migration of the vapor to the reactor walls, (ill) a parametric study 
using the boundary layer code to optimize the performance characteristics of 
the Westinghouse reactor, ' iv ) calculations relating to the collection effi- 
ciency of the new AeroChem reactor, and (v) final testing of the modified LAPP 
code for use as a method of predicting Si(i) droplet sizes in these reactors. 
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I. INTRODUCTION 


The goal of this program Is to develop mathematical models and computer 
codes that will predict product distribution In chemical reactors designed to 
convert gaseous silicon compounds to condensed-phase silicon. Emphasis is 
placed on collecting silicon as a liquid on the reactor walls and other col- 
lection surfaces. The two reactor systems being studied are the Westinghouse 
SiCl*/Na reactor 1 in which Si(£) is collected on the flow tube reactor walls, 
and the new AeroChem reactor in which Si(t) droplets formed by the SiCl,./Na 
reaction are collected by a jet impingement method. 2 

Conceptually, it is useful to visualize a flow reactor of the Westinghouse 
type 1 as having two distinct parts: 

(1) An "upstream" section in which gaseous silicon halide (e.g., SiCl k ), 
with argon and hydrogen as diluents, entrains into a jet of sodium vapor, 
mixes rapidly due to turbulence, and reacts to yield the required silicon- 
containing product gas stream. 

(ii) A "downstream" section, with cooled reactor walls ( *■ 1700 K) , through 
which the hot ( 3500 K) developing flow of product gases passes. Separation 

(and collection) of silicon occurs in this section by the passage of condensed 
silicon droplets and uncondensed vapor, through the turbulent boundary layers 
onto the reactor walls. 

The situation is depicted in Fig. 1. Because of the difficult nature 
of both problems, our analysis also has two distinct parts. 

The first part of the analytical work, using LAPP (AeroChem' s rocket plume 
code), considers the problem of the open, turbulent, axisymmetric jet with fi- 
nite rate, multiple-step chemical reactions. Recent modifications of LAPP have 
extended its capabilities to i..«-lude a detailed kinetic description of silicon 
droplet formation via homogeneous nucleation and hence can be used to compute 
the streamwise histories of droplet size and number concentration along the 
reactor.* On the other hand, the silicon separation/collection problem in the 
"downstream" section is dominated by the presence of walls and is better 


Test runs are currently being made and the modified code is nearly debugged. 
The solutions from LAPP will be an important input in the solution of the 
"downstream" problem. 
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FIGURE 1. SCHEMATIC OF THE VARIOUS REGIMES WITHIN A SILICON REACTOR 


handled with the well-known GENMIX (Version 4A) computer program 3 for the 
solution of boundary layer type (i.e., parabolic) problems. This latter 
problem is the major subjec. of discussion in this report. 

II. WALL DEPOSITION OF SILICON 

The case of silicon vapo r deposition, caused by the various convective- 
diffusion processes in the developing, turbulent flow field encountered over 
most of the "downstream" reactor section, is analyzed here using operating 
conditions typical of the current Westinghouse reactor design. 1 Existing 
Westinghouse analyses 1 of the silicon separation/collection problem use the 
following rather simple model. By considering a differential axial segment 
of the flow as a control volume, these analyses express mass, momentum, and 
energy balances in terms of the net efflux of these quantities across the 
control volume boundaries. The details of the boundary layer structure near 
the reactor walls are ignored by the utilization of empirical heat and mass 
transfer coefficients (Nusselt numbers), appropriate to a given developing 
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turbulent pipe flow, in particular, the streetwise variation ; the Nusselt 
number (for given Reynolds and Prandtl numbers) is taken from the experimental 
results of Boelter et al, 4 pertaining to the case of air flowing through a 
pipe with a right-angle bend entrance. Furthermore, the effects of condensa- 
tion on the wall heat and mass transfer rates are accounted for via a cor- 
rection factor based on a fully-developed (i.e., no streamvise variations in 
flow properties) turbulent pipe flow analysis involving silicon vapor (no 
droplets). Thus condensation is treated solely as an augmenter of heat and 
miss fluxes. Using this one-dimensional, constant property analysis it was 
concluded that the silicon vapor collection efficiency of the proposed 8 m 
long reactor would be about 801. 

Thus, while the Westinghouse analysis outlined above serves as a useful 
tool for preliminary reactor design, a unified analysis (such as the one we 
adopt) will ultimately be needed to make reliable design predictions. A prop- 
er consideration of th A dominant deposition mechanisms and the structure of 
the developing turbulent boundary layer is essential before measures can be 
taken to improve the silicon collection efficiency of a given reactor. In 
what follows we discuss exactly such a model for the case of silicon vapor 
deposition. Undoubtedly, the presence of condensing silicon droplets in the 
developing turbulent flow field will alter the silicon collection efficiency 
since the turbulent, Soret, and Brownian diffusional behavior of particles 
can differ greatly from that of vap molecules. 5 This latter aspect of the 
problem will be the subject of future reports. 6 ** 8 Finally, it might be worth 
mentioning that when the details of the silicon deposition behavior are 
adequately modeled and computerized (as we have done) they will enable the 
designer to produce an optimum reactor design with minimum experimental 
effort. 


III. SILICON VAPOR DEPOSITION/COLLECTION MODEL 

We start with the assumption that all chemical reactions have reached 
completion, and hence the product gas stream is chemically inert, just prior 
to entering the downstream section (region II, Fig. 1) of the reactor. Thus, 
in this section the model formulated will describe the behavior of the various 
convection and diffusion processes responsible for depositing silicon on the 
reactor walls. The physical problem involves an analysis of the hot, turbulent. 
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developing flow field within e circular pipe with cooled walla. Due to thu 
large temperature variation between the reactor core flow and the walls, it 
is essential to Include the effects of non-uniform fluid properties and the 
possible formation of silicon droplets due to condensation (caused by cooling 
silicon vapor). However, for simplicity, in the present analysis we restrict 
attention to the case with no condensation. That is, silicon is assumed to 
exist only in the vapor state. 

Yet another simplification is rendered in this analysis by viewing the 
overall product gas mixture as an effective binary system, with silicon as 
one species and all other gases (mainly NaCI, Ar, H a ) as the only other species 
(referred to subsequently as the "carrier" gas). The molecular weight and 
other properties of the carrier gas are taken to be a weighted average of the 
properties of the individual components. On this basis, the diffusional be- 
havior of silicon vapor through the surrounding gas can be uniquely described 
in terms of the appropriate mass diffusiviti.es (or mass transport coefficients). 9 
On the other hand, the thermal diffusivity (or heat transport coefficient) of 
the overall product gas mixture will be nearly equal to that of the carrier gas 
alone, since silicon concentrations are typically low. In the case of silicon 
vapor deposition, the relative importance of heat and mass transfer processes 
inside a typical reactor is nearly equal. That is, the Lewis number for the 
system is close to unity. The effect of turbulence, as treated in the model 
considered, lies mainly in enhancing the diffusive transport of heat, mass, 
and momentum.* The enhancement due to the so-called turbulent counterparts 
of the abovementioned dif fusivities will be assumed to be such chat the 
"effective Lewis number" remains close to unity. That is, turbulence enhances 
heat and mass transfer rates equally. However, according to current belief 
and experimental data the enhancement of momentum transfer rates due to tur- 
bulence is somewhat less (a fact represented by setting the turbulent Prandtl 
number'*’ equal to 0.9). 


The diffusive transport of momentum (or vorticity) is analogously controlled 
by the kinematic viscosity (v), or momentum diffusivity. 

The Prandtl number is the ratio of momentum diffusivity to thermal 
diffusivity. 
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Ic is clssr, therefore, thst a good description of turbulent transport 
coefficients is pivotal to the success of the present model. For now, a mix- 
ing length approach to turbulence modeling has been adopted. The chief draw- 
back of this approach is that it relies heavily on an empirical determination 
of several constants, some of which might vary from one flow situation to an- 
other. This handicap will have to be overcome by the use of more sophisticated 
(differential) models of turbulence. This latter class of turbulence models, 
which describe the mixing length via differential equations (rather than alge- 
braic ones, as done here), can, of course, be blended into the present formula- 
tion. However, we believe that since turbulent heat/mass transfer predictions 
using even the most elaborate approaches sometimes incur errors of up to 15Z,* 
retaining the mixing length model "'.aid at least give us the advantage of great 
er simplicity with relatively little loss of achievable accuracy. Moreover, 
the mixing length model has been extensively studied and tested in pipe flow 
configurations against both experimental data 3 and the predictions of more 
sophisticated models.’ 0,11 

In view of the abovementioned assumptions and recognizing that the develop 
ing pipe flow processes of interest are characterized by comparable changes in 
the radial and axial directions (i.e., two-dimensional, axisymmetric) one may 
adopt the Patankar and Spalding 3 turbulent boundary layer formulation for the 
conservation of overall mass, momentum, energy (or stagnation enthalpy) and 
silicon vapor species.* No attempt will be made to reproduce these well- 
established parabolic, partial differential equations here since they have 
been fully detailed in Refs. 3 and 12. Instead here we choose to discuss 
aspects not immediately obvious from their discussion. 

A, THE PATANKAR-SPALDING FORMULATION 

Implicit in the governing equations used by these authors 3,11 are the 
following assumptions: 


It can be estimated from the magnitudes of the various dif fusivities that 
the velocity, temperature, and species boundary layers will all develop at 
nearly the same rate. The "fully-developed" state is not reached until 
some 40-t0 diam downstream from the pipe inlet. 
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(1) The density fluctuations caused by turbulence are of Insignificant 
importance compared to the flux contributions of other fluctuating quantities 
(e.g. , velocity, temperature, concentration). 

(ii) The dissipation rate of turbulent kinetic energy (due to velocity 
fluctuations interacting with the mean velocity) is negligible compared to 
the viscous dissipation rate of mean kinetic energy. 

Under these assumptions, it is possible to reduce the otherwise complicated 
turbulent boundary layer equations to a form similar to the better-understood 
steady, compressible, laminar boundary layer equations, provided one defines 
suitable "effective" values for the momentum, thermal, and mass diffusivities 
discussed earlier. These have been defined by the following flux-gradient re- 
lationships (assumir^, Newtonian behavior) : 


T 


U 


3u 

eff 3y 


(shear 3tress) 


= -r c 

" h,eff p 3y 


(heat flux) 


3Y, 

1 " = -r — - 

J i i, ef f 3y 


(mass flux of species i) 


where U e j£ is the effective viscosity, the ratio of effective thermal 

conductivity (A ) to the constant pressure specific heat of the carrier gas 
erf 

(C ), and T, the product of the mean gas density (p) and the effective 
P i » f 

mass diffusivity (D^ e ^) of silicon vapor in the carrier gas. u, T, and Y^ 
denote the time-averaged streamwise velocity, absolute temperature, and mass 
fraction of silicon, respectively. 

Each of the effective transport coefficients above is the sum of a 
laminar (or molecular) contribution and a pseudo, turbulence-induced con- 
tribute n. The latter may be greater by some two orders of magnitude in 
fully turbulent (i.e., nearly inviscid) regions, such as the reactor core 
flow. However, closer to the reactor walls turbulence energy is rapidly de- 
pleted due to the dissipating influence of molecular viscosity, and transport 
is affected mainly by the laminar (molecular) mechanism. As a result the 
turbulent parts of the transport coefficients can vary significantly across 
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the reactor cross-section. The motivation for introducing a "mixing-length," 
or any other turbulence model, is to describe this variation realistically. 

In the present analysis the following radial variation of mixing length, 
from the reactor axis (or centerline) to the wall, was prescribed: 


£ 


j 16 • y>f 

\ Ky [ 1-exp (-y + /A + )], 0 ^ y ^ ~ 


(co-i layer region) 


(near wall region) 


where y is the distance from the wall, 6 is the thickness of the turbulent 

+t 

boundary layer, and A, K, and A are taken to have the constant values 0.09, 
0.435, and 26.0, respectively. Physically, such a distribution of mixing 
length distinguishes regions of high local Reynolds number, y = yu*/v, (or 
nearly inviscid flow) from those regions near the wall which are governed by 
the progressively increasing effect of molecular viscosity (i.e., a "damping" 
of the mixing length with decreasing y + , according to Van Driest's 14 exponen- 
tial law) . Note that this permits a viscous sublayer region next to the reac- 
tor walls so that the influence of even weak turbulent fluctuations is felt. 
Just outside the viscous sublayer the mixing length obeys an undamped, "defect 
law" behavior (i.e., £ * Ky) . 

Using the above mixing length distribution it is possible to express the 
effective viscosity as: 


eff 


u turbulent + U laminar 




t ii + 

It has been shown ' more recently that A is a sensitive function of pressure 
gradient aid wall blowing or suction. This parameter determines the thick- 
ness of the "viscous sublayer" egion. 

^ u^ is the so-called "friction velocity" defined here in terms of the local 
shear stress by the relation u^ = (x/p ) 1 / 2 
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where |-r — 1 is the magnitude of the time-averaged velocity gradient normal to 
dy 

the reactor wall. Since this quantity decreases away from the wall (becoming 
zero at the boundary layer edge) turbulence is operative only in regions in 
which the velocity gradient exceeds a small critical value. If, in any flow 
region the turbulence velocity scale |£ falls below a certain fraction of 
the local velocity, this quantity is set to the lcc^l prevailing velocity. 
Such a circumstance tends to occur in the fully-turbulent core flow region 
of the reactor. 


B. COMMENTS ON THE NUMERICAL SOLUTION PROCEDURE 

Since the governing equations to be solved are parabolic partial differ- 
ential equations, their solution requires the specification of appropriate 
initial and boundary conditions. In the present study the initial conditions 
at the pipe entrance were specified as uniform (corresponding to centerline 
values) profiles of u, T, and Y across the pipe cross-section. The domain 
of integration was taken to be the region between the pipe centerline and 
the wall. Boundary conditions on u, T, and were specified at the wall as: 

u » 0 

T = T (constant) - 1700 K 

w 

Y. = Be~ Ay/T w =» 0* 

l 

The last condition is a consequence of assuming that silicon vapor is in 
equilibrium at the wall (i.e., the partial pressure equals the saturated 
vapor pressure) . The centerline boundary conditions on these variables werp 
taken to correspond to the symmetry conditions: 


3u 

9y 




0 


* 


A and B are constants with the values 46710.0 K and 7.3166 x 10 10 
respectively. 


N m 


-2 


s 


Using the above initial and boundary conditions, GENMIX produced solutions 
at specific downstream stations along the reactor by marching forward in steps 
(whose size was proportional to the local boundary layer thickness) . At each 
step the solution to the coupled system of conservation equations is carried 
out using a fully implicit, 6-point, finite difference algorithm (based on an 
integral approach that ensures conservation, rather than the usual Taylor se- 
ries expansion). The algorithm uses a substitution method for solving the re- 
sulting system of coupled algebraic equations (with a tri-diagonal matrix), 
thus avoiding the time-consuming and unreliable operation of inverting the 
coefficient matrix. It was found in this study, as in others, that GENMIX is 
extremely efficient as a code. For instance, the full solution for a 9 m 
reactor length with some 80 to 100 radial grid points is completed in about 
17 s (CP time on CDC 7600 computer). 

The remarkable efficiency of GENMIX is really a result of the many dif- 
ferent time-saving features built into the code. While the program chooses 
its owr forward marching step size, the cross-stream grid spacing is specified 
as input. In the present problem, the grid points were so chosen as to be 
unequally spaced over the pipe radius, with a greater resolution capability 
near the wall, where gradients were steep. Furthermore, it must be noted that 
the solution in GENMIX is not carried out in physical space (x,y) but rather 
in the transformed von Mises 1 ® coordinates (x,w) . In this latter coordinate 
system, the normalized stream function o> always provides a fixed integration 
domain: 0 ^ w ^ 1. In general, the advantage of working in this transformed 

space is the greater solution accuracy and ease that results from imposing 
boundary conditions at fixed extremities. However, in the solution of con- 
stant radii pipe flows (such as the present case), this advantage is not fully 
realized since the two extremities remain fixed even in physical space. Of 
course, an alternative solution tactic may have been to solve the problem be- 
tween the wall and the edge of the developing boundary layer. In this latter 
case, working in the transformed stream function ordinates may be a definite 
advantage. 

A major disadvantage of the GENMIX solution procedure may be its handling 
of the near-wall flow. The reactor wall is a point of singularity in the 
(x,u) space because the flow velocity there is zero and it multiplies the 
highest derivatives of the transformed equations. This disconcerting feature 
also tends n- ender the transformation process back to physical space somewhat 


9 


Inaccurate for grid points close to the wall. In order to circumvent some of 
the difficulties associated with this singlar point, the Patankar-Spalding pro- 
cedure generates solutions for the region next to the wall using a turbulent 
Couette flow analysis (further simplified using a constant property assumption) . 
Thus one does not obtain a rigorous finite difference solution of the two- 
dimensional boundary layer equations in the region close to the wall. However, 
this device greatly enhances the numerical solution efficiency since wider grid 
spacings can be used. Moreover, by arranging the region of Couette flow very 
close to the wall, the error due to this approximation may be rendered acceptable. 

C. RESULTS AND DISCUSSION 

Of the numerous results that were generated by GENMIX we focus attention 
here on those which especially pertain to the mechanism of silicon vapor deposi- 
tion and/or help us improve the performance of a given silicon reactor. 

1. Cross-Stream Profiles 

The structure of the developing, turbulent flow in the "downstream" 
section is characterized by the velocity, temperature, and silicon mass frac- 
tion profiles across the reactor tube radius. Figure 2 shows the variations 
in normalized (with respect to the centerline values) velocity, temperature, 
and concentration vs the normalized radial distance from the pipe centerline. 

The operating conditions and pipe geometry were chosen to match the Westing- 
house reactor design. 1 As expected from the prescribed mixing length dis- 
tribution discussed earlier, it is possible to discern at least three regions 
of distinctly different behavior. Near the pipe centerline (r/R = 0) the 
fully turbulent core flow tends to even out all gradients due to rapid mixing. 

On the other hand, near the wall (r/R = 1) the steepest gradients are estab- 
lished under the influence of the molecular dif fusivities for momentum, heat, 
and mass competing with residual (damped) turbulent dif fusivities. Within 
this "viscous sublayer" the turbulent fluctuations decay rapidly (according 
to the Van Driest hypothesis) as the wall is approached. Thus, in the imme- 
diate vicinity of the wall, the final deposition of silicon may occur pri- 
marily by molecular mechanisms (i.e., a nearly laminar flow condition). 

Joining these two extremes is the outer turbulent (defect law) layer where 
transport occurs mainly due to diffusive turbulent fluctuations, characterized 
by a mixing length that varies linearly with distance from the wall. Note 
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that in this region, unlike the viscous sublayer, there is negligible damping 
due to viscosity. Rather, the length and velocity scales of turbulenc ' are 
smaller than in the core flow. 

Although in reality the abovementioned zones blend continuously into 
each other, it has been customary to vlev turbulent pipe flows in terms of at 
least three different layers (with respect to the velocity profile) 1 *: the 

laminar sublayer, the buffer layer, and the fully turbulent layer. Such a 
model is especially true of fully developed pipe flows. However, even Fig. 2 
(for a developing flow) reveals these zones. The profiles are nearly linear 
next to the wall, implying that the laminar sublayer is a region of constant 
fluxes of momentum, heat, and mass. It is also interesting to note that the 
thicknesses of these three laminar sublayers are different, with the tempera- 
ture and mass fraction layer thickness being equal and larger than the velocity 
layer thickness. Such a behavior is dictated by the input values of the laminar 
Prandtl and Schmidt numbers. In this case Pr ■ Sc * 0.7 was used. On the 
other hand, the buffer layer region is also seen to be one of nearly constant 
fluxes, even though the fluxes here are less than those in the laminar sub- 
layer. Based on the prevailing wall shear stress, these two zones are expected 
to correspond to regions of the velocity profile that lie in the ranges 0 ^ y + 

^ 5 (laminar sublayer) and 5*^y + ^ 30 (buffer layer). Both these zones fall 
within the viscous sublayer. Elsewhere in .he boundary layer, the fluxes de- 
crease with distance away from the wall and the final boundary layer thick- 
nesses are nearly all equal. This is a consequence of the assumption that 
the effective Lewis and Prandtl numbers are nearly unity. 

2 . Streamwise Variation of Fluxes 

Figure 3 shows the variation of the non-dimensional momentum, heat, 

and mass fluxes along the reactor length, X (non-dimensionalized with respect 

to the reactor diameter, D). The flow and other conditions are identical to 

those used ^or Fig. 2. The following definitions of the Stanton numbers for 

heat (St.) and mass (St ) transfer, and the skin-friction coefficient (C c ), 
h m r 

were employed : 


St 


h 
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, u b (h b 
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suggests that under these conditions a Reynolds analogy of the type: 

St - C f /2 

might be valid within the reactor since the effective Prandtl (and Schmidt) 
numbers are close to unity and the pressure variation in the downstream direc- 
tion is small due to a low inlet velocity (in this case U . , . » 20 m s~ l ). 

c, inlet 

In the other cases studied, it was found that greater deviations from this 
Reynolds analogy occurred with increasing inlet velocity due to more signifi- 
cant pressure gradients. Table I summarizes some of these results. It 
should be pointed out that "bulk” values of properties differ from their 
corresponding centerline values (denoted by the subscript "c"). This is 
clear from the variation of l^/U^ with downstream distance. Typically, the 
centerline values monotonically exceed the bulk values in the region of de- 
veloping flow and become constant only when "fully developed" flow conditions 
are established.* 


It has been shown 9 that in a pipe, as the developing shear layers merge, 
exceeds (U c ) fully developed P rior to settling down to a state of equality. 
Such details of shear layer interaction, however, cannot be predicted by the 
present mixing length model. 
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/U, is tl>e ratio of centerline to bulk velocity. 
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3. Silicon Collection Efficiency 

The ratio of the mass of silicon deposited on the reactor walls, due 

to the boundary layer convective-diffusion processes, to the total mass of 

silicon entering a given reactor may be defined as the "collection efficiency." 

This ratio can readily be calculated by integrating the downstream variation 

of the silicon mass flux to the walls, j", over the reactor length of interest. 

w 

Figure 4 shows the results of such a computation. It is seen that the silicon 
collection efficiency, for a given inlet flow velocity, increases as the reac- 
tor length increases. Initially, this increase is slower than in the aft re- 
gions (as expected from the mass flux variation shown in Fig. 3). Furthermore, 
Fig. 4 reveals the important effect of changing flow residence times on the 
reactor's performance. As expected, faster through-flows reduce the collection 
efficiency drastically since the flow time is much less than the time for 
effective silicon vapor diffusion to the reactor walls. Charts of this type 
can help the designer produce an "optimum reactor," at least from the stand- 
point of silicon separation and collection. 



FIGURE 4. COMPUTED VARIATION OF SILICON VAPOR COLLECTION 
EFFICIENCY OF A REACTOR WITH LENGTH AND INLET VELOCITY 


*» 
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IV, SILICON COLLECTION VIA JET IMPINGEMENT 

The new AeroChem SiCl„/Na reactor 2 will employ a considerably different 
Si collection scheme from that used by Westinghouse. 1 In this new reactor a 
Jet impingement method of collecting Si(£) droplets will be tested. It is 
expected that this collection method will be effective in collecting the very 
small Si(£) droplets produced via the SiCl„/Na reaction and in effectively 
separating the Si (A) from the NaCl produced by the reaction. The reactor and 
Si collection crucible to be tested shortly are shewn in Fig. 5. Na(g) and 
SiCl*(g) react in a graphite chamber, 8 cm in diam and 16 cm long. The reac- 
tion temperature can be adjusted by heating the reactor, and initial tests 
will be done with ** 1700 K walls; the gas within the reactor is xpected to 
be at S 2300 K. Flow rates of reactants will be adjusted to maintain a pres- 
sure of 0.5 atm (360 Torr). At the bottom of the reactor a sir pie converging 
nozzle with a throat diameter of 1.2 cm will allow the reacted NaCl(g)/Si(£) 
droplet mixture to exhaust into a large vacuum vessel held at a pressure (for 
initial testing) of **10 Torr. (Since the products of reaction are condensi- 
ble and only a very small amount of Ar is used to keep the SiCl* delivery line 
shielded from the hot product gases only a very small pumping capacity is re- 
quired after the initial pumpdown. ) The jet of product gas/droplets will im- 
pinge on an Si(A) surface at the heated crucible. This will result in a nor- 
mal shock standing above the Si(£) surface and in a post-shock layer at high 
temperature and moderate pressure through which the droplets must pass to 
reach the surface. Using standard isentropic flow and normal shock tables and 
the operating conditions given above, the values of pressure, temperature, den- 
sity, and gas velocity in the reactor, the jet and the post-shock layer given 
in Table II are obtained. 

Calculations have been done to find the minimum Si droplet size which 
might be expected to penetrate the post-shock layer and impinge on the sur- 
face, i.e., be collected by impingement. The first ttep in the calculation 
is to find the post-shock layer thickness. To do this one assumes that the 
thickness is uniform and determined by the fact that the radial mass flow be- 
neath the shock must be equal to the mass flow of gas in the jet. One thus 

has, for the p' 5t-shock thickness, l , , 

sh 

l • 1/2 p u r /(u p ) 

sh J j j s s 
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TABLE II 

TENTATIVE OPERATING CONDITIONS FOR THE AEROCHEM REACTOR 



Reactor 

Jet 

Post-Shock Laver 

Pressure, Torr 

360 

9.5 

110 

Temperature, K 

2300 

810 

2180 

Density, gm ml - * 

1.37 (-4) 

1.00C-5) 

4.0(-5) 

Velocity, cm s“* 

0 

1.36(5) 

3.3(4) 

Mach No. 

0 

3.02 

0.47 

where p is density. 

u is velocity. 

and the subscripts j 

and s refer to the 


jet and post-shock regions. The radius of the jet, r ^ , is 2.07 times the 
exit radius of the nozzle, or 2.5 cm. Using the values for density and 
velocities, of Table II, we find 


£ , = 1.3 cm 

sh 

The calculation of stopping distances for droplets entering the post- 
shock layer poses a major problem. Values for stopping distances calculated 
two ways are given in Table III. A glance at this table reveals that the 
Knudsen numbers, N^, (the ratio of molecular mean free path to droplet 
radius) for the very small droplets of interest (diam ■ 0.5 pm) are very 
large. This raises a problem since the usual Stokes relationship for com- 
puting stopping, even with corrections for molecular slip, will not 

be trustworthy for such large values of N,, . This standard relationship 

Ml 

is the formula based on Stokes flow: 


P d u 

5 «= - P - P -J 

st 18 p 


x SCF 


where p g is the viscosity of the post-shock gas, p g = 1.1 * 10 -3 g s~ l cnT l , 
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TABLE III 

STOPPING DISTANCES FOR DROPLETS ENTERING THE POST-SHOCK LAYER 
Layer thickness = 1.3 cm 


Droplet 
Diam, pm 

Knudsen 

No. 

Stokes 
Stopping 
Distance, cm 

Stopping Distance, 
based on Molecular 
Collisions, cm 

0.01 

1220 

0.04 

0.03 

0.05 

240 

0.19 

0.58 

0.10 

120 

0.39 

1.57 

0.50 

24 

2.0 

12.4 


p and d are the droplet density (2.4 g ml l ) and diameter, respectively, and 
P P 

the Stokes-Cunningham slip factor is 

SCF - 1 + 1.26 N + 0.4 N v expC-l.l/N^ ) 

Table III includes values for S gt obtained using this formula. Physically, 
the trouble with such a formula is the assumption that a particle, as it 
moves through the gas, carries with it an enveloping cloud (boundary layer) 
of gas with it. For >>10 this will not be the case and a better approach 
would be one in which particle kinetic energy losses result from individual 
particle/gas molecule collisions assuming that the gas molecules in front of 
the oncoming particle are not disturbed by the particle's approach (i.e., a 
"free-molecule" approach rather than a "continuum" approach). 

If is the kinetic energy of a particle of mass m^, then the loss of 
kinetic energy suffered by collision with a single gas molecule of mass m 

S 

(considered on the average to be stationary) will be (with averaging perform- 
ed over the surface of the particle) 

AE - -2 m E /m 

P g P P 
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The number of such collisions suffered per cm of travel will be (ird n )/4 

P 8 

where n is the number of gas molecules per ml. The collisional stopping 
8 

distance S ^ is then given by 


S coll 




X 


2 m 

g 

Trd a n m 
P 8 P 


where E *■(«■ 1/2 m u 2 ) is the particle kinetic energy as it enters the post- 
P P ] 

shock layer and E^ is the final particle kinetic energy. The values of E^ 
specified cannot be zero and for Table III have been taken to be 


E * = 1.0 x 10“ k T 

P s 

i.e., 6700 times the mean thermal translational energy of gas molecules at the 

post-shock temperature. Fortunately, the calculation does not depend strongly 

on the choice of E^ ; a change of a factor of 10 in its value changes the values 

of S by a factor of 2.3, (i.e.. In 10). 
coii 

From Table III it is seen that S , , > l . for d =0.1 um. Even if the 

coll sh p 

Stokes stopping distance S is picked as a conservative value it is found 

that droplets with d = 0.3 um will reach the liquid surface by impingement. 

P 

Thus it is found that the jet impingement method is capable of collecting 
very small droplets and, if sue’ droplets are produced. Si collection effi- 
ciencies will be high. The sizes of droplets expected will shortly be examined 
using the modified LAPP code now in its final testing stages and operating 
conditions will be sought which will yield drop . . .s with d^ ~ 0.1 um. 

V. CONCLUSIONS 


A model and computer code have been developed to describe the behavior of 
the turbulent, developing boundary layer flow in the "downstream" section of 
a tubular reactor, where the crucial silicon separation/collection is affected 
The present analysis necessitated modification of the GENMIX-4A code 3 in order 
to achieve solutions to the coupled system of parabolic, partial differential 
equations for conservation of overall mass, momentum, energy, and species. The 
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turbulence model adopted was of the mixing-length type, with wall damping 
according to a modified Van Driest hypothesis. In particular, the following 
objectives have been achieved: 

(i) Silicon vapor deposition processes can now be analyzed with due 
regard to the structure of the developing turbulent flow that 
prevails over most of a reactor's length. 

(ii) The modified GENMIX code now provides a computational capability 
that should be useful to the designer in optimizing a given 
reactor, with a minimum of experimentation. 

(iii) A reliable basis for assessing the silicon vapor collection 
efficiency of a given reactor has been established. 

Finally, it must be mentioned that Soret transport of silicon vapor was 
intentionally excluded from this analysis, since it is expected to be small 
in the absence of condensation. The effect, however, will become important 
when silicon droplets are considered in future analyses . 4-8 ’ 17 

Calculations have also been performed which indicate that the "jet 
impingement" Si collection method to be used with the new AeroChem SiCl<,/Na 
reactor is capable of collecting very small Si droplets. 

VI. PLANS 


During the next quarter, boundary layer calculations which will include 
Si condensation and thermophoretic droplet deposition 3 on the reactor walls 
will be made to more accurately determine Si collection efficiencies for the 
Westinghouse reactor, under a variety of operating conditions. The modified 
LAPP code will be used to predict Si droplet size distributions in both the 
Westinghouse and AeroChem reactors. 

VII , NEW TECHNOLOGY 

No reportable items of new technology have been identified. 
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